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ABSTRACT 

In this paper we present evidence that the H ll regions of high luminosity in disk 
galaxies may be density bounded, so that a significant fraction of the ionizing photons 
emitted by their exciting OB stars escape from the regions. 

The key piece of evidence is the presence, in the Ha luminosity functions (LPs) of 
the populations of H ll regions, of glitches, local sharp peaks at an apparently invariant 
luminosity, defined as the Stromgren luminosity (Lstr), = Lstr = lO^*^'^ (± lO"'^) 
erg s~^ (no other peaks are found in any of the LFs) accompanying a steepening of 
slope for Liia> Lstr- This behavior is readily explicable via a physical model whose 
basic premises are: (a) the transition at Lhq- = L^tr marks a change from essentially 
ionization bounding at low luminosities to density bounding at higher values, (b) for this 
to occur the law relating stellar mass in massive star-forming clouds to the mass of the 
placental cloud must be such that the ionizing photon flux produced within the cloud 
is a function which rises more steeply than the mass of the cloud. Supporting evidence 
for the hypothesis of this transition is also presented: measurements of the central 

1/3 

surface brightnesses of H ll regions for Lna < i^str are proportional to L^^^ expected 
for ionization bounding, but show a sharp trend to a steeper dependence for Lua > -^str, 
and the observed relation between the internal turbulence velocity parameter, sigma, 
and the luminosity, L, at high luminosities, can be well explained if these regions are 
density bounded. 

If confirmed, the density-bounding hypothesis would have a number of interesting 
implications. It would imply that the density-bounded regions were the main sources 
of the photons which ionize the diffuse gas in disk galaxies. Our estimates, based on 
the hypothesis, indicate that these regions emit sufficient Lyman continuum not only 
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to ionize the diffuse medium, but to cause a typical spiral to emit significant ionizing 
flux into the intergalactic medium. The low scatter observed in Lstr, less than 0.1 mag 
rms in the still quite small sample measured to date, is an invitation to widen the data 
base, and to calibrate against primary standards, with the aim of obtaining a precise, 
~ 10^ Lq, widely distributed standard candle. 

Subject headings: ISM: H ll regions - galaxies: individual (NGC 157, NGC 3631, 
NGC 6764, NGC 6951, NGC 6814, NGC 5194, NGC 7479, NGC 3359) - galaxies: 
spiral 
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1. Introduction 

The basic theory for modehng a gaseous re- 
gion round a hot star was first given by Zanstra 
(1931) for planetary nebulae, and applied to H II 
regions by Stromgren (1939), who quantified the 
relation between the radius of the ionized zone 
and the temperature-luminosity of the central 
exciting star. He showed that, in a uniform 
medium, the transition layer between fully ion- 
ized and neutral gas will be thin compared to 
the radius, a structure since termed a Stromgren 
sphere. This will occur in a large placental neu- 
tral cloud, of sufficient dimension to absorb all 
the ionizing photons (those in the Lyman con- 
tinuum in the case of atomic hydrogen) , in which 
case the H II region is "ionization bounded" . The 
case where the cloud is not big enough to absorb 
all the ionizing radiation was first dealt with by 
Hummer k, Seaton (1964) for planetary nebulae, 
which are then termed optically thin. An H II 
region formed in these circumstances, where the 
cloud radius is less than that of the Stromgren 
sphere, is termed "density bounded" or "matter 
bounded" . 

There is ample evidence that H il regions are 
not homogeneous in density. Measurements of lo- 
cal electron densities via line ratios are typically 
two orders of magnitude higher than mean elec- 
tron densities estimated using diametral emission 
measure (for a clear recent example see Rozas, 
Knapen, & Beckman 1996), explained if a re- 
gion comprises knots of high density embedded 
in a lower density plasma. The fractional vol- 
ume dense enough to contribute measurably in 
emission lines such as Ha has been termed the 
"filling factor" [e.g., Osterbrock 1989). An H II 
region with this structure will form a Stromgren 
sphere provided that the local density variations 
are on scales which are small compared with the 
H II region diameter, and that the mean density 
varies little on the scale of the diameter. This 
condition appears to hold well for regions over 
a wide range of Ha luminosities, up to a criti- 



cal value, which we will be able to quantify from 
observations, as described below. The regions in 
this range are ionization bounded. At higher lu- 
minosities our evidence, which will be described 
in this paper, points towards an increasing ten- 
dency for the photon output from the ionizing 
stars to overflow the cloud in which they have 
formed. 

One type of evidence comes from the measure- 
ments of a change in slope of the luminosity func- 
tion in Ha of complete populations of H II regions 
in the set of nearby spirals. Previous detections 
of this change have been reported in the litera- 
ture, notably in earlier work by Hodge, and by 
Kennicutt and co-workers (Hodge 1987; Kenni- 
cutt 1984; Kennicutt, Edgar, & Hodge 1989); 
and most recently by McCall and co-workers 
(McCall, Straker, & Uomoto 1996; Kingsburgh & 
McCall 1998) as well as ourselves (Rozas, Beck- 
man, & Knapen 1996). It is interesting to note, 
however, that the change in slope varies little 
from object to object. We will elaborate this 
point during the development of the present pa- 
per. 

Statistical studies of the relation between the 
Ha luminosities, Lho;, of complete samples of H II 
regions in nearby large spirals, and their volumes, 
(Cepa & Beckman 1989, 1990; Knapen et al. 
1993; Rozas, Beckman, k. Knapen 1996), show 
that to a first approximation these are propor- 
tional. Taken complete description of the 
observations this would imply, at least statisti- 
cally, two properties: (1) The H II regions obey 
a single physical regime, which should be that of 
ionization bounding, and (2) the densities of the 
clouds in which the H II regions are formed differ 
little from cloud to cloud within a galaxy, and 
from galaxy to galaxy. We will see that at very 
high luminosities departure from this first order 
behavior is dominant. Illustration of these points 
may be found in a number of articles in the liter- 
ature, where the Ha luminosities of individually 
measured H II regions have been plotted, loga- 
rithmically, vs. the cubes of their radii (Cepa & 
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Beckman 1989, 1990; Knapcn et al. 1993; Rozas 
Knapen, & Beckman 1996; Rozas et al. 1999). In 
these articles, it is shown that for eight galaxies 
there is a close to linear relation between the Ha 
luminosity and the volume of the H II regions in 
the range 37.5 < log Lhq < 39.5, as pointed out 
in Cepa &: Beckman (1990). This implies that 
the product of the filling factor, the electron den- 
sity and the ionized hydrogen density does not 
change very mTich over this range, although in 
all the galaxies there is a clear tendency for it to 
increase at the highest luminosities. 

It is not entirely ruled out that the filling fac- 
tor and mean density might vary independently 
with a relation which leaves their product invari- 
ant, but this would amount to a "conspiracy", 
and it is more reasonable to assume that none of 
the cited parameters varies strongly. The general 
picture is of H ii regions as spongy structures, 
which may indeed have fractal characteristics as 
suggested by Elmegreen (1997), but with a de- 
gree of clumping, and an average density which 
do not depend strongly on luminosity. However 
in all the plots of log Lhq vs. cited above, 
there is a significant trend to steeper gradients 
at higher luminosities, which means that here the 
densities and/or the filling factors of the regions 
must be increasing. 

The evidence which we present below points 
to a rather sharp change in the properties of the 
H II regions occurring at an Ha luminosity which 
appears to vary very little from galaxy to galaxy. 
We will show that one would in fact expect, on 
reasonable physical grounds, the most luminous 
regions to be producing more ionizing photons 
than can be absorbed in their placental clouds, 
whereas this should not be true globally for the 
less luminous regions. We can predict, from a 
model in which there is a transition between ion- 
ization bounding at lower luminosities, and den- 
sity bounding at higher luminosities, that there 
should be an accumulation of H II regions around 
the luminosity of the transition (which we have 
termed Lstr, ^-c, the Stromgren transition), be- 



tween the two regimes, and this accumulation 
is in fact found "glitch" in the Ha lumi- 

nosity function of the galaxies observed. It is 
by no means easy to reproduce these observa- 
tions with alternative hypotheses. Two aspects 
of the density-bounding phenomenon would be 
of particular interest. Since the most luminous 
regions are those from which we predict the high- 
est proportion of Lyman continuum (Lyc) pho- 
tons are escaping, these regions are candidates 
for the principal ionizing sources of the diffuse 
inster stellar medium in large spirals, and may 
contribute significantly to the intcrgalactic ioniz- 
ing field. Also the transition between ionization 
bounding and density bounding appears to take 
place over a narrow range in Lrq, and while the 
sample of results presented here is not large, this 
point is worth following up as a high-luminosity 
feature of this kind would have its uses as a pow- 
erful, precise, and non-transient standard candle. 
However the detailed interpretation of the sur- 
face brightness data is not obvious and further, 
spectroscopic, observations will be important to 
follow up the arguments presented here. 

In Sect. 2 we present the evidence about 
the change in properties of the H II regions, 
in Sect. 3 we present scaling relations which 
support the model of density bounding for the 
highest-luminosity regions, and a set of alter- 
native models aimed at explaining the observed 
"glitches" in the luminosity functions of H II re- 
gions in the galaxies observed. In Sect. 4 we dis- 
cuss the implications of our models for the diffuse 
Ha emission from the warm ionized medium. In 
Sect. 5 we summarize our conclusions, outline 
the observations required to deepen our under- 
standing of the transition phenomenon, and sug- 
gest applications to the correction of the global 
star formation rates (SFRs) in galaxies, and to 
the measurement of intcrgalactic distances. 

2. The evidence for a change in physical 
regime 
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Fig. 1. — Three-dimensional surfaces: best fits to the observed data representing the numbers, N, 
diameters, d, and Ha luminosities, -Lna, of the H II regions in six galaxies. The projections in the log 
N-log Liia plane are luminosity functions, LFs with effective binning of 0.1 dex. The broad peaJiS and 
declines to low Lho are caused by observational selection: low-luminosity regions are increasingly difficult 
to measure reliably below log Lhci ~ 37, though all samples are complete to below log Lho = 38. The 
sharp peaks at Lhq; ~ 38.6 are presumably due to the transition from ionization bounding to density 
bounding (see text for details). 



5 



51 




Fig. 1. — Continued 

2.1. The luminosity functions 

In a set of published papers stretching back 
over the past decade we have pubhshed luminos- 
ity functions (LFs) in Ha for the H II regions 
in spiral galaxies, the majority of which were in 
absolutely calibrated fluxes (Cepa &: Beckman 
1989, 1990; Knapen et al. 1993; Rozas, Beck- 
man, & Knapen 1996). It was while inspecting 
these LFs, with the addition of the LF from Rand 
(1992), that we found the feature which led us to 
postulate the transition from ionization bound- 
ing to density bounding. All the LFs showed 
an abrupt change in slope accompanied by a 
"glitch" : a jump upwards in the value of the LF, 
at a luminosity close to log Lh^ = 38.6 {Lua in 
ergs~^). In Fig. 1 we show an improved analysis 
of the LFs originally given in Rozas, Knapen, & 
Beckman 1996, for the galaxies NGC 157, NGC 
3631, NGC 6764, and NGC 6951, in Knapen et 
al. (1993) for NGC 6814, and in Rand (1992) 
for M51, in which the principal change has been 
an effective rebinning into bins of 0.1 dex in lu- 
minosity. The new treatment of the data uses a 
three- dimensional surface fit of log Lna against 
number count (log A'^) and diameter, which in 
Fig. 1 has been projected into the log Lna-log 
N plane with a resolution of 0.1 dex. The data 
show an absence of peaks in each LF except for 
the peak at the transition luminosity which, as 



NGC 5951 




seen, is accompanied by a change in the slope of 
the LF to steeper values at luminosities higher 
than that of the transition. 

We carried out a specific check on NGC 157 
to see how the glitch peak is affected if instead 
of using a limiting surface brightness contour to 
define the integrated brightnes of each H li re- 
gion we applied a cut-off at a constant fraction 
of the peak luminosity. This technique ("PPF") 
was described in Kinsburgh Sz McCall (1996); it 
avoids problems of S:N ratios at the boundaries 
of the H II regions, and should be superior to the 
limiting surface brightness method ("FTP") for 
sets of geometrically homologous regions. If there 
is a trend to varying internal brightness gradient 
with luminosity, the PPF technique is liable to 
imprecision. In the present context the most lu- 
minous regions have very bright centres and ex- 
tended haloes, so that PPF tends to underesti- 
mate their luminosities. This was confirmed with 
the LF of NGC 157 which fell off more rapidly 
at high Lffa with PPF than using our normal 
technique, i.e. FTP. In spite of this, the glitch 
showed up, at the same luminosity but at a sta- 
tistically less significant level, due to the smaller 
numbers of H II regions involved. 

To check that a peak is plausibly due to a 
truly collective property of the H ii regions, and 
is not merely a statistical fluctuation, we took 
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the mean value, N(m.), of the number of H II re- 
gions in the two bins to either side of the peak, 
subtracted this from the value N(max) of the 
number of regions at the peak, and computed 
{N{max) — N{m))/N{max)^^'^ . Values between 
1 and 2 were found for each of the galaxies, with 
a mean of 1.7. Thus each individual peak is just 
significant, but if we were looking at a single 
galaxy and had no other kinds of observational 
evidence we would not, perhaps, want to draw 
any striking conclusions from the presence of a 
peak no higher than 2a. However the peaks oc- 
cur in a number of galaxies, and the chance that 
this occurs as the only such feature, at closely 
similar luminosity, for purely statistical reasons, 
is very small, less than 1 part in 10^, so it would 
be very hard to dismiss the peaks as statistical 
artifacts. The fact that the peaks and gradient 
changes occur at the same luminosity within 0.07 
mag, for galaxies whose distances cover a range of 
over a factor three, also argues strongly against 
an observational effect (e.g., limitation in angular 
resolution) as their cause. 

A qualitative explanation for the presence of 
the glitches in the LFs in terms of the transition 
from ionization bounding to density bounding is 
as follows. In a hypothetical set of star- forming 
clouds with identical cloud masses but with star 
clusters of different integrated ionizing luminosi- 
ties forming inside them, some of the clusters 
may be massive enough to emit more than suffi- 
cient ionizing photons than are required to fully 
ionize their placental clouds. In that case, what- 
ever the intrinsic LF of the resulting H ii regions 
(measured in ionizing photons) might be, the LF 
observed in Ha would show a spike at the upper 
bound of the observed luminosity range, because 
any regions with extra ionizing photons would 
lose the excess by density bounding. In this ex- 
treme, and not of course realistic, case we would 
not find a glitch, but a delta function in the LF. 
Cloud masses are not, however, identical. In or- 
der to yield a glitch two entirely plausible physi- 
cal conditions are required: one is that the mass 



of massive stars rises, statistically, with the cloud 
mass; the other is that this rise is sufficiently 
steep that at a critical cloud mass, and above 
it, there are sufficient Lyc photons to ionize the 
whole cloud. Wc will see in Sect. 3 that the ev- 
idence which docs exist in the literature is fully 
consistent with these two conditions. 

2.1.1. Modeling of the LFs 

Here we show, using a set of fairly schematic 
models, that the process of density bounding 
ought to give rise to the form of the LF observed 
in Ha, and that the other hypothetical causes 
modeled for the observed behavior are in no way 
as plausible. The model starts by generating an 
artificial LF using a population of 1000 H ll re- 
gions, distributed in bins of uniform logarithmic 
width; we chose 0.2 dex for the illustrations here 
as this is our preferred plotting scale for the ob- 
servations. The number of regions per bin, n, is 
given by 



dN 



dlog L 



n{L 



Ha J 



Ha 



-^Ha 
A' ' 



(1) 



where Lho is the luminosity in Ha, and the 
index ^ can be varied, but in our figures we will 
keep this at an observationally typical value of 
-1.5. The constant A' is adjusted in order to 
normalize the number of regions at log L^ict — 
38 erg s~^ in the 0.2 dex bin to 100. The form 
of equation (1) is the most direct way to repre- 
sent the luminosity function, but conventionally 
either differential or logarithmic forms are used 
for convenience. These are respectively: 



and 



diV(LHa) 

dN 



L 



Ha 



A 



dLfia 



T 1—a 
-^Ha 

A' ' 



(2) 
(3) 



C^log ^Ha 

where a = 1 — ^ and A = A' /\oge. The form 
(3) will be used in the development of Section 3. 

With the basic function of equation (1), we 
can then apply processes which simulate the dif- 
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Fig. 2. — Computed models showing how the 
phenomenon of density bounding above a criti- 
cal luminosity can lead to the observed glitch plus 
a change of slope in the luminosity functions of 
the H II regions of a disk galaxy (see text for 
details). Each contour plot corresponds to an in- 
trinsic slope of the initial LF (values are shown). 
The dashed plot show an LF with slope -1.5 prior 
to modeling the effect of density bounding. 

ferent physical effects expected to affect the LF. 
In the first place we show a simulation of the ef- 
fects of density bounding. The model assumes 
that the observed luminosity of a region in Ha is 
reduced to Lna-, above the critical value of Lho 
(Lstr) at which an H II region is just fully ion- 
ized, by a factor which depends on the ratio of 
the true ionizing luminosity (Lj) to the ionizing 
luminosity which produces Lhq = -Z^str = -^i,Str- 
This is a fair parametric representation of what 
is predicted when the ionizing flux overflows the 
H II region. The formula for this is just: 

^Ha = X (-^)^ (4) 

-i^i.Str 

where A; < and L^^ is the H II region lu- 
minosity in Ha if all the ionizing photons are 
absorbed. Given the dependence of the observed 
Ha luminosity of a region on its intrinsic ioniz- 
ing luminosity, we can transform numerically to 
obtain the corresponding LF in Ha. Applying 



this to the part of the LF, given by equation (1) 
in the range Lh^ > Lstr (for < i^str, L-^a 
= L^^), with a smoothing function of Gaussian 
type with cr = 0.1 dex in luminosity (to corre- 
spond to the binning of the observational data) 
we obtain curves of the type shown in Figure 2. 
We can see here that the general functional form 
of our observations is reproduced, notably the 
glitch in the LF as well as the increase in (nega- 
tive) slope for Lfja > ^str- We can also see that 
the luminosity of Lstr is almost independent of 
the slope, ^, of the initial trial function, which 
represents the intrinsic slope of the LF and re- 
flects essentially the IMF slope of the massive 
exciting stars within the region, i.e., the glitch 
peak luminosity should depend little on the ini- 
tial mass function (IMF) if the break is due to 
density bounding as we suggest. 

To examine another possible cause of the break 
and the glitch in the LF we assumed that the H ii 
regions contain an admixture of dust, causing ex- 
tinction in Ha. Dust is a well-explored feature 
of the insterstellar medium, and the main doubt 
is to what extent it survives within the extreme 
conditions inside an ultra-luminous H II region. 
For the purpose of this exercise we assumed a 
constant dust-to-gas ratio, and uniform mixing, 
which is the case most favorable to a dust ex- 
tinction hypothesis as the cause of the LF break. 
This was treated via a modification of the intrin- 
sic Ha luminosity L'^^ given by eq. (1) (n(L^^) 
= L'yi^/A), to a resultant observed luminosity 
L-Ra-i according to the expression: 

Lna = L'hq (5) 

where r is an optical depth in Ha given by 

T = B r, V = 4/37rr"^ and the constant B takes 
values which normalize the extinction as follows: 
at a radius of r = 100 pc a region has an intrinsic 
luminosity L^^ = 10^^ erg (chosen as a rep- 
resentative value from the observations), which 
is proportional to its volume, V (again a depen- 
dence taken from the observations, to first order). 
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and B is then varied at will to give values of r in 
simple decimal intervals. 



10000.0 F 




og (erg/s) 

Fig. 3. — Computed models of H II region LFs 
with constant gas-to-dust ratio. Each continuous 
curve corresponds to a value of integrated dust 
optical depth (at the Ha wavelength) for a region 
with emitted Ha luminosity of log Lhq = 38 (erg 
s~^). The dotted curve is the unextinguished LF. 

In Figure 3 we can see the results of this mod- 
eling exercise for different values of the r at log 
Lua = 38 erg s~^. The results show qualitatively 
very different LFs from what we observe. If the 
optical depth is small the effect is to curtail the 
LF but not to alter its slope. If the optical depth 
is larger, the effect is to curtail the LF with a high 
peak and a drop to zero at a cut-off luminosity, 
without altering the slope of the main LF. 

Although the result of this exercise with dust 
suggests that dust extinction is unlikely to be 
the cause of the LF we observe, it does point to 
a possibility which at least qualitatively might 
produce an LF not unlike that observed. If dust 
is mixed within the H II regions up to some limit- 
ing luminosity, but then in those of higher lumi- 
nosity is expelled, it may be possible to modify 
the curves in Fig. 2 in such a way as to displace 
some of the H ll regions accumulated in the rising 
tail of the LFs in Fig. 3 to higher luminosities, 
yielding an LF of the form observed. While this 
scenario is not ruled out physically, it would give 



rise to a dependence of the central Ha surface 
brightness of the H ll regions on radius which is 
very different from that observed, as we will ex- 
plain below, and this is sufficient reason for us 
not to accept it as the explanation for the LF 
observations. 

Finally we modeled numerically another pos- 
sible scenario for the break in the LF: could it 
be caused by the overlapping of H II regions of a 
certain size and above? In constructing the H II 
rgion catalogue, one of the difficulties that there 
is to overcome is that many H II regions appear 
to overlap on the image. We adopted the solution 
proposed in Rand (1992) of counting each peak 
in Ha as representing a single H ii region. The 
flux of each H ii region was then estimated by 
integrating over the pixels which could be rea- 
sonably attributed to a given region. One will 
undoubtedly miss a number of H ii regions that 
are too weak to be detected in the vicinity of 
stronger emitters close by. This will influence 
the lower end of the LF but is not a significant 
factor in the determination of the true LF at the 
higher luminosity end (Rand 1992). Anyway, to 
check directly if the overlapping cloud causes the 
break in the LF, we applied a new model. First 
we assumed the function in eq. (1) as the intrinsic 
LF. To take care of the weaker regions in a way 
corresponding closely to reality, we used eq. (1) 
between 10^° erg s~-^ and 10^''' erg s^^ and below 
this limit we used a flat distribution function (^ 
= in eq. (1)) following the observed distribu- 
tion in M31 by Walterbos & Braun (1992). In 
the computations we allowed the luminosity of 
a region merged with others to reach an upper 
limit of 10*^^ erg s~^, corresponding to the lu- 
minosity of the regions formed by the mergers. 
The distribution was discretized into 300 bins of 
equal logarithmic width, and the volume of H ll 
regions per bin was calculated by simple summa- 
tion, normalizing the system so that an H II re- 
gion of radius 100 pc had an Ha luminosity of log 
Lua = 38 (erg s^^) and calculating the volumes 
in the first instance at infinite spacing between 
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regions, i.e., with no overlaps. The luminosity 
per bin was proportional to the summed volume 
of regions in the bin, as we found observationally 
for the regions below log Lua = 38 (erg s~^). The 
total volume of all the regions is V^. 

In the model wc put the regions whose inte- 
grated volume is Vc into a enclosure of volume 
V, so the overall filling factor is Vc/V. The ini- 
tial distribution is considered in finite luminosity 
(i.e., volume) bins, indexed as i, j, where i and 
j can both vary from 1 to 300. Interaction of the 
set of clouds, one in bin i and the other in bin j, 
is first considered. 




0.1 

36 37 38 39 40 

og Lho (erg/s) 

Fig. 4. — Models of the LF in Ha in which in- 
creasing degrees of overlapping and absorption 
take place between H ii regions (see text for de- 
tails). Although we place little stress on the 
quantitative results, qualitatively we can see that 
the effect of assimilation of regions within others 
is to reduce the LF slope at high Lua, and not 
to increase the slope as observed. 

The probability per unit volume, Pj, of a small 
cloud absorption within a larger, is the filling fac- 
tor of the clouds in bin j, which is the same as 
using the rule that when two clouds are sited 
with the distance between their centers less than 
the radius of the larger, the smaller is absorbed 
into the larger, and the resulting H ii region has 
the sum of the volumes (and luminosities) of the 
two. The mean initial spacing is determined by 



the chosen value of Vc/V. We obtain for Pj: 



where Nj is the initial number of clouds within 
j, so the number SNf of clouds affected is just: 

5Nf = NiNjVj. (7) 

We then adjust the luminosities in proportion 
to the revised numbers of clouds in each bin via: 

Ni' = Ni- 5Nf 

N/ = Nj - 6Nf (8) 
Nk' = Nk + SNf 

in which the total volume and luminosity are 
conserved, while the regions affected by merging 
are transferred to the appropriate bin of higher 
luminosity. This procedure is carried out in 
order of increasing bin luminosity, so that the 
merged regions participate in further mergers if 
they are close enough to other regions. In a sin- 
gle pass combining all pairs of bins the new LF 
is computed, modeling the mergers in this way 
we assume that H ii regions are not optically 
thick, so that only genuine three-dimensional 
overlaps, not projected overlaps on the plane of 
the sky, need to be taken into account. As will 
be deducible considering Fig. 4, assuming opti- 
cal thickness would only accentuate the disagree- 
ment between this type of models and the ob- 
served LFs. In Fig. 4, it is clear that not only 
is there no evidence of a glitch, but the negative 
slopes of the IMFs of the merged models become 
less steep at high luminosities, in direct contrast 
to the observations. While it might be possible to 
find some type of clustering model which would 
yield an LF similar to those observed, its inputs 
would have to be very finely tuned, and would 
be unlikely to reflect the realistic astrophysical 
conditions. 

In the light of the comprehensive theoretical 
study of H II region LFs, notably in the Galaxy, 
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but with commentary on other galaxies, by Mc- 
Kee & WiUiams (1997, hereafter MW) we should 
add here that the effect of density bounding on 
the LF in Ha which we claim here must function 
independently of almost any detailed consider- 
ation of the IMF of the stars which ionize the 
regions, or its temporal evolution. The trunca- 
tion of the IMF at high stellar mass, adduced 
by MW to explain an observed sharp steepening 
of the LF slope of H II regions does not conflict 
with the results presented here for two reasons. 
Firstly the phenomenon explained by MW occurs 
at an Ha himinosity almost an order of magni- 
tude below the glitch we find. In the present 
study, this would be very difficult to distinguish 
from the rather rapid change in slope just above 
log i^Ha= 37, due to the increasing incomplete- 
ness of our statistics at low luminosity, but in 
any case it is clear that one is dealing with two 
different real observable effects. Secondly a trun- 
cation of this type does not readily give rise to a 
glitch in the LF, nor can it yield as a natural con- 
sequence the observed sharp increase in internal 
brightness gradients and central surface bright- 
ness with luminosity, described below in Section 
2.2. Thus neither the observations nor the model 
of MW relate to the same phenomenon as the one 
we describe here. 

2.2. Central surface brightness and inter- 
nal brightness gradient 

Further evidence for a change in physical regime 
at the Ha luminosity Lstr comes from measure- 
ments of the surface brightnesses of H II regions. 
The observations used are the same as those for 
the LFs: photometrically calibrated CCD im- 
ages of galaxies taken through an appropriately 
redshifted Ha filter, with stellar continuum sub- 
traction via a neighboring off-band filter. All 
observations were taken with the 4.2-m William 
Herschel Telescope on La Palma, and details of 
the observations and the reduction procedure can 
be found in Rozas, Beckman, & Knapen (1996). 
The absolute Ha fluxes and mean radii were 



found for several hundred H ii regions per galaxy. 
The largest, most luminous, and isolated were se- 
lected for profile measurements, as these subtend 
angles large enough to obtain reliable profile fits; 
a minimum radius of 5 pixels was chosen. This 
corresponds to an Ha luminosity log Lua = 38 
in the furthest galaxy observed. 

The numbers of regions per galaxy which could 
be used for these measurements was 30-35. Two 
characteristic parameters were derived per re- 
gion: the central surface brightness and the mean 
internal surface-brightness gradient (this latter 
was also treated in Rozas, Castafieda, & Beck- 
man 1998). 

In Fig. 5a we show central surface brightness 
1/3 

as a function of L^^^, the cube root of the 
Ha luminosity. For spherical ionization-bounded 
regimes with constant density and filling factor, 

Sc should be proportional to the radius of a re- 
gion, and as luminosity will be proportional to 

1/3 

volume, Sc should vary linearly with L^^. We 
can see from Fig. 5a that for regions with lu- 
minosities below Lho ~ Lstv this relation holds 
rather well, but for luminosities above this value 
there is a clear and sharp departure from propor- 
tionality. In Fig. 5b we show a brightness gradi- 
ent parameter, obtained by dividing the central 
surface brightness of a region by its linear radius 
projected in the plane of the galaxy, plotted ver- 
sus Liia- We can see that below Lna = -^str this 
parameter varies very little, within the limits of 
error, while above this luminosity it tends to rise. 
Although the scatter also increases, the trend to 
higher values is clear, and the ratio of this gradi- 
ent parameter above Lstr to the constant value 
below Lstr reaches 10 in the steepest cases, and 
is greater than 3 in most regions with log Lhc* > 
39. The change in behavior seen in Figs. 5a and 
5b could not be due to the finite resolution of the 
images. It is easy to show that for a hypothetical 
set of regions with constant intrinsic brightness 
gradient the observed parameter, if degraded by 
instrumental resolution, would tend asymptoti- 
cally to a constant value at large radii, i.e., at 
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high luminosity. This is just the converse of what 
we measure: a tendency to a constant value at 
low luminosity. If on the other hand the intrinsic 
gradient were in fact decreasing with increasing 
luminosity it would take a "conspiracy" to flatten 
the plot at low Lhq (i.e., an inadequate deconvo- 
lution folded with an intrinsic fall would need to 
cancel neatly, yielding a level plot), and the re- 
sult would disagree even more strongly with the 
observations at high LjHa- It would be very sur- 
prising to find the observed change of gradient 
at a similar luminosity for all the galaxies sam- 
pled if it were due to a resolution effect, since 
they are at quite different distances, so the lin- 
ear resolution on the galaxy disk corresponding 
to the angular resolution of the telescope differs 
from object to object. The observed changes in 
the value of the gradient above Lhq = ^Str are 
way outside the most pessimistic error bars, so 
that Figs. 5a and 5b offer prima facie evidence 
of a change in physical parameters close to Lua 
= Lstr- The regions below this luminosity obey 
a linear volume-luminosity relation, while those 
above it in general do not. This break in the ge- 
ometrical properties of the regions cannot be ex- 
plained plausibly in terms of the collective prop- 
erties of the ionizing stars. It is easy to show, 
and it was shown by Stromgren himself, that for 
an H II region forming within a medium of con- 
stant uniform density (and, we should add here, 
constant filling factor) the volume ionized is pro- 
portional to the rate of emission of ionizing pho- 
tons, so that the geometrical distribution of the 
brightness within a region will not depend on the 
number of its ionizing stars, provided that they 
are concentrated in a cluster near the center (cer- 
tainly true for all the regions studied in this sec- 
tion, which have radii of well over 100 pc). The 
observed changes must be due to the behavior of 
the surrounding gas. We will discuss in Sect. 3 
of the paper how these changes might occur for 
the observed regions. 

A further point of interest here is that the con- 
stancy of the brightness gradient parameter, and 



the dependence of the central surface brightness 

1/3 

on Ljj^ below Lhq = i^str are not consistent with 
the dust model which we described in Sect. 2.1 
as a testable alternative scenario for the observed 
behavior of the LF. Dust extinction increasing 
linearly with the radius of the H II region, and 
causing a notable bump in the LF close to log 
Lua = 38.6, would have an effective value of opti- 
cal depth THa (at Ha) close to 0.6 at our normal- 
ization value of log Lhq = 38, while thq would 
be 0.4 at log Lua = 37.5, and 0.9 at log Lua 
= 38.5. The surface-brightness gradient would 
show an easily detectable systematic fall between 
37.5 and 38.5, by a factor of order 2, rather than 
remaining constant as observed. If above Lua = 
Lstr the observed increase in the surface bright- 
ness were due to dust being dispersed in these re- 
gions (perhaps by extra radiation pressure), the 
gradient would not increase by factors of 3 to 10 
from its value at LHa= -^Str, but at most by a fac- 
tor 2, and in any case it would never rise above 
the asymptotic value at small Lhq, a rise which 
is dramatically observed. Very similar consider- 
ations hold for the central surface brightnesses. 
We will see in Sect. 4 that the diffuse Ha ob- 
served in the disks of our galaxies offers more 
direct evidence that ionizing radiation is in fact 
leaking out of the luminous H ii regions, and this 
tells against major dust extinction in these re- 
gions. 

2.3. Internal velocity dispersion of the 
H II regions 

The evidence discussed in the present sub- 
section does not refer to the same set of galaxies 
as those treated in the bulk of the paper, because 
the kind of data required, Fabry-Perot mapping 
of a galaxy in an emission line, is in no way as 
easy to obtain as the narrow-band imaging on 
which the majority of our deductions are based. 
We have, to date, analyzed Fabry-Perot maps in 
Ha of the disks of only two galaxies, obtained 
in a program of disk-wide kinematics, not aimed 
principally at the properties of H II regions. How- 
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ever the statistical results on the internal veloc- 
ity dispersions of the H II regions observed are 
of such relevance to the present discussion that 
it is important to include them here. They were 
made with the TAURUS instrument on the 4.2-m 
William Herschel Telescope on La Palma. Details 
of the observations and reduction are presented 
elsewhere (Rozas, Sabalisck, Sz Beckman 1998). 
Wc present the data here for two galaxies: MlOO 
and MlOl. We had absolute flux calibrations for 
the regions in MlOO, obtained by our group in a 
parallel context (Knapen 1998), but for MlOl we 
have no such absolute fluxes, nor are they avail- 
able in the literature, so we give the results in 
relative flux units. This still allows us to infer 
the essential relations between flux, Lhq, and in- 
ternal velocity dispersion, a, for the regions of 
MlOl. 

Figure 6a is a plot of log Lh^ ^s. log a for the 
principal components of the brightest H ii regions 
in MlOO. Only the brightest 10 %, some 200 of 
the almost 2000 regions whose absolute fluxes 
and diameters had been determined (Knapen 
1998), were bright enough to yield velocity pro- 
files of high enough signal to noise ratio to give 
well-measured values of a from single-frame ex- 
posures. We also have a lower limit cut-off in 
a of 10 km s~^, due to the combined uncertain- 
ties in the instrumental profile, and the natural 
plus thermal broadening contributions. At first 
glance Fig. 6a is little more than a scatter dia- 
gram, and any attempt to put a linear, or even a 
polynomial or spline fit through the points yields 
little of value. One might be tempted to con- 
clude that the horizontal spread is a little less 
than the vertical spread, so that a line of high 
gradient, perhaps akin to the prediction based 
on the assumed virial equilibrium of the H ll re- 
gions (Terlevich & Melnick 1981), 

log Lua = 4 log a + c, (9) 

could be brought to fit the points, but there is 
no real justification for such a fit. However there 
is one locus of regularity in these points, namely 



the upper envelope of minimum a, which can be 
well fitted by a straight line of slope 2.6 (±0.14). 
It is worth comparing this with the linear fit to 
the log Lhq vs. log a plot for the H II regions 
in MlOO by Arsenault, Roy, & Boulesteix (1990) 
using a similar observational technique to ours. 
In their case they selected the regions of great- 
est surface brightness, rather than an envelope 
of minimum a. This coincidence is very signifi- 
cant physically, as wc will show in Section 3. In 
Fig. 6b we show a similar log Lna-log cr plot for 
MlOl; although the luminosity is not on an ab- 
solute scale, as explained above, wc can sec here 
also the linear envelope to what would otherwise 
be a scatter diagram. The best fit measured slope 
to this is 2.55 (±0.15), and the display of points 
in the graph is similar to that in Fig 6a, with a 
slight trend for the a values to occupy a wider 
range at lower luminosities. 

Our interpretation of the observations in Figs. 6a 
and 6b is that, in general, the internal velocity 
dispersions of the H II regions as shown in their 
Ha emission line profiles are not in virial equilib- 
rium. Whatever the mechanisms for transferring 
energy to the ionized gas (and it is beyond the 
scope of the present article to consider these) it 
is reasonable to assume that the emission line 
width, for an arbitrary region, will be enhanced 
by the presence in the gas undamped eflFects of 
the outflowing winds from massive stars, or even 
of supernova explosions. In the region NGC 
604, in nearby M33, Sabalisck et al. (1995) 
showed the morphological complexity of the ve- 
locity field. The principal velocity component 
for a more distant region, observed via our tech- 
nique, is a luminosity-weighted sum, which must 
incorporate any discrete burst of input energy, 
and should not, in general, fit a virial equilib- 
rium relation. In this scenario, regions showing 
the lowest value of a for a given luminosity will 
be those whose emission profiles are least affected 
by non-virialized energy injection, so that regions 
closest to virial equilibrium should lie on, or near, 
the low-(T envelope of the log Lna-log o relation. 
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Physically, it is plausible that the most massive 
regions will damp most rapidly the bursts of in- 
jected energy, so that in the full log Lna-log a 
diagram an increasing fraction of regions will be 
found on the virial envelope with increasing mass 
(and hence luminosity). The linear envelopes in 
Figs. 6a and 6b can thus be identified as the 
loci of those regions whose Lna-c relation should 
obey the conditions of the virial theorem. 

It is therefore very significant that the slopes 
of these envelopes are ~ 2.5 and not 4 as pre- 
dicted for virialized systems (see, for example, 
Terlevich & Melnick 1981; Tenorio-Tagle, Munoz- 
Tunon, & Cox 1993), since the value of 4 fol- 
lows from assuming the conventional picture of 
ionization-bounded regions. The regions on the 
envelope in Fig. 6a all have luminosities greater 
than Lgtr, so we would expect them to be den- 
sity bounded, and we will show below that a slope 
close to 2.5 is in the range predicted for virialized 
regions which are density bounded. Arsenault, 
Roy, &; Boulesteix (1990) performed a similar ex- 
ercise to ours for the H II regions of MlOO, but 
selected a small number of regions of maximum 
surface brightness for their log Lna-log cr plot. 
The slope of their graph, 2.6, confirms our sce- 
nario, as we have shown above that the regions of 
highest surface brightness are found in the range 
-^Ha > Lstr, where the regions are most massive, 
where the virial relation should hold in the ion- 
ized gas, and where our theory predicts density 
bounding. 

In this section we have presented evidence for 
a change in the physical properties of the H II re- 
gions in spirals in the very high luminosity range 
(which were termed "giant H ii regions" by Ken- 
nicutt, Edgar, & Hodge 1989), above Lho = 
Lstr = 10^^'^ erg s~^ in Ha. We have argued, 
based firmly on observation, that this may mark 
a "phase change" from ionization bounding to 
density bounding. In the next section we argue 
from scaling calculations that the occurrence of 
this type of change is physically plausible. 



3. Why should very luminous H II regions 
be density bounded? 

3.1. The luminosity function 

In Sect. 2.1 we outlined what would be re- 
quired of the relation between the mass in high- 
mass stars within a star-forming cloud and the 
mass of the cloud itself so that the H ii regions 
with luminosities above a specified value were 
density bounded, while those below this value 
were ionization bounded. The relation is that the 
mass in ionizing stars should, with cloud mass, 
grow quickly enough for the ionizing flux to rise 
more quickly than the cloud mass. In the very 
schematic model we present below we assume, 
to simplify the argument, that the average cloud 
density does not vary, an assumption which holds 
well, observationally, for regions up to the tran- 
sition. The model will not fail physically if this 
condition is not well maintained to higher lumi- 
nosities, but the simplicity of the quantitative 
inferences will be lost. 

Let the ionizing flux, Li, from an H II region 
depend on the total stellar mass, M*, in ionizing 
stars within the region via 

Li = k{M^Y, (10) 

where A; is a constant. We can set observa- 
tional constraints on a. From the semi-empirical 
study by Vacca, Garmany, &: Shull (1996) we find 
that the rate of emission of Lyman continuum 
photons from OB stars rises approximately as the 
square of the stellar mass. The index falls from 
values higher than 2 at early B and late O to val- 
ues close to 2 at 03. For a uniform "Salpeter" 
IMF slope, and no physical cut-off in stellar mass 
at the upper mass end with increasing placen- 
tal cloud mass, the weighted mean index for a 
young star cluster would be a little over 2, and 
this would be the appropriate value for a. If, 
however, there were a fixed physical stellar mass 
limit, rather than a statistical limit for each clus- 
ter, a would be unity or close to unity. Although 
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evidence on this point is not abundant, Massey 
& Hunter (1998) measured with HST the pho- 
tometric properties of the most massive stars in 
the cluster R136 which ionizes the very luminous 
H II region 30 Doradus in the LMC. They found, 
inter alia, that the IMF slope appears to retain 
its Salpeter value of -2.35 to the highest masses, 
and that there is no obvious evidence for an up- 
per mass cut-off; they find stellar masses ranging 
well over 120 M©, through 130 Mq and up to 
150 M0. Without generalizing unduly, or be- 
ing especially optimistic, it is reasonable to infer 
that a, though not necessarily greater than 2, 
will generally be higher than 1, and wc will say 
provisionally that a should be between 1 and 2. 

Now we let the mass in stars contribut- 
ing significantly to the ionization depend on the 
placental cloud mass Mci according to 

M^ = jM^„ (11) 

where j is a constant. There is not very 
much direct evidence on the value of e. Lar- 
son (1982) reported, from a compilation of ob- 
servations within the Galaxy, that the highest 
mass of a young cluster varies as the placental 
cloud mass to the power 0.43. If we take the 
IMF slope as -2.35, the mass of stars contribut- 
ing to the ionizing photon flux rises with an in- 
dex of ~ 0.43x2.35; we derive this value of 2.35 
as a sum of 4.35 from the mass-weighted inte- 
gral of the Salpeter function, and -2 which takes 
into account the effective contribution of mass to 
ionizing luminosity. Using Larson's result thus 
yields a value for epsilon of 1.01, effectively of 
unity. 

The condition that the rate of production of 
ionizing photons fully ionizes, and overflows the 
placental cloud of a young cluster for clouds 
whose masses (and therefore H il region lumi- 
nosities) exceed a specific value is found by com- 
bining (10) and (11), and is just 

ae > 1. (12) 



The empirically inferred values for a and e de- 
rived above readily satisfy this condition. Al- 
though, as we have seen, direct evidence related 
to a and e is sparse, there are two indirect ways 
to estimate ae from the observations presented 
in the present paper, both of which confirm that 
a e must be greater than unity, one using the 
gradient of the Ha LF, and the other from the 
slope of the virial envelope of the log Lna-log cr 
diagram. 

The first estimate comes from the ratio of the 
slopes of the H ii region LF in Hq above and 
below the critical luminosity Lna = i^str- 
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Fig. 7. — Diagram showing how the integrated 
ionizing flux from the luminous density bounded 
H II regions in NGC 157 can be estimated by 
subtracting the observed LF for log > 38.6 
from the extrapolated LF obtained by extending 
the observed function below log Lhq = 38.6 to 
higher values. This estimate indicates that the 
density-bounded regions leak sufficient Lyc pho- 
tons, in this galaxy, to ionize the diffuse medium 
without the need for other sources of ionization. 

We can express the LF for Lhq < -^Str as: 
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dNi{LHa)=AL-^dLua (13) 
and for Lho; > ^Str as: 

dN2{LHa) = BL-ldLna, (14) 

where A and B arc constants. Inspection of 
Fig. 7 shows that this parameterization is a good 
representation of the data. In the range Lho 
< Lstr (below the transition), the functional ex- 
pression for the ionizing flux is the same that for 
the Ha flux in equation (13) (as these H ii re- 
gions are ionization bounded); so that from (10) 
we can write 

dNi{M^) = bM-^dM^, (15) 
where 6 is a constant, and 

5=l + a(/3-l). (16) 
Combining (15) with (11) yields: 

dNi{Mci) = cM^'dMcu (17) 
where c is a constant, and 

77 = 1 + 6(5-1). (18) 
Substituting (18) and (16) into (17) we have: 

diVi(Mci) = cM-(^+'°(^-^))dMci. (19) 

Equation (19) is the mass function of the 

clouds below the transition (Ljia < ^str)- 

If clouds with luminosities Lh^ > ^str are 
fully ionized, the slope of their mass function is 
also that of their luminosity function. This slope 
is the index 7 of cq. (14). At Lho = -^Str the 
slopes of the mass functions below and above Lstr 
should be equal, as there is no reason to predict 
a discontinuity in the cloud mass function. Then 
we can write 



7 = l + ea(/3-l), (20) 

which is 

ae =(7 -!)/(/? -1). (21) 

Prom eq. (21) we deduce that the condition 
ae > 1 implies 7 > /3, i.e., that the slope of the 
LF in the range Lna > -^Str should be steeper 
than the slope in the range Lho < Lstr- We 
can see in Fig. 1 that the LFs for all the galax- 
ies observed do meet this condition, i.e., for each 
galaxy 7 is greater than /?. In Table 1 we give 
numerical values for /? and 7 for eight galaxies 
taken from our own work and from the litera- 
ture. Since e ~ 1, we could use eq. (21), and the 
measured LF in Ha for a galaxy to obtain an em- 
pirical value for a. Another way to derive a is by 
weighting the data in Vacca, Garmany, & Shull 
(1996) by the stellar IMF. Equating the two val- 
ues enables us to derive the IMF slope. Using 
this technique, the values found lie in the range 
2 to 2.5 for the LFs of the galaxies measured in 
the present paper. However the method is not 
intrinsically very accurate, as it requires the ra- 
tio of two obscrvablcs, each of form (a:; — 1) where 
X ^ I. Another, relatively more reliable way to 
obtain the IMF slope uses the gradient of the log 
-^^Ha~log cr envelope in the velocity-dispersion di- 
agram, as we will now explain in Section 3.2 

3.2. The internal velocity 

The virialization of the gas motions within 
H II regions has been proposed (Terlevich & Mel- 
nick 1981) not so much as a working hypothesis, 
rather desirable case, which would allow us 
to measure the gravitational mass of an emit- 
ting region. We showed in Sect. 2.3 that the 
majority of the H II regions in the two galax- 
ies whose log Lfja^log cr diagrams are plotted 
are not in virial equilibrium. The Ha flux they 
emit is a combination of individual sources within 
each region, some of which have non-random mo- 
tions with relatively high velocities. Although 
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NGC 157^1^ 


38.66 


1.31 


2.18 


NGC 3359(2) 


38.60 


1.44 


2.10 


NGC 3631(1) 


38.54 


1.51 


2.16 


NGC 5194(3) 


38.63 


1.52 


2.39 


NGC 6764(1) 


38.60 


1.30 


2.11 


NGC 6814(4) 


38.69 


1.67 


3.39 


NGC 6951(1) 


38.49 


1.33 


2.74 


NGC 7479(5) (disk) 


38.65 


1.25 


-2.10 



Table 1: The measured Ha luminosities, Lstr erg , at the peaks of the glitches, for the seven galaxies 
where data of high quality yielding complete LFs in Ha are available. Luminosities were measured using 
standard stellar calibrators, the distance to the objects were taken on the basis of a value of 75 km s~i 
Mpc~i for iifo, assuming that a measured radial velocity is due exclusively to the Hubble flow In the 
succeding columns are the slopes of the LFs in Ha below and above Lstr, /3, and 7 respectively. The 
mean measured uncertainty in /? is ±0.05 and in 7, ±0.15. (^^ Rozas, Beckman & Knapen 1996; Rozas, 
Zurita & Beckman 2000; Rand 1992; Knapen et al. 1993; Rozas et al. 1999 



the mass fraction in these, typically expansive, 
motions may not be high, a luminosity-weighted 
profile incorporates them, and will always tend to 
have a super-virial value for its half-width. We 
argued in Sect. 2.3 that the lower envelope in a 
to the log Lfja-log a distribution should corre- 
spond to the virialized regions, and showed that 
the regions found on this lower bound in a (which 
is also an upper bound in Lna), will in general 
be density bounded. One corollary of this is that 
the predicted power relationship between lu- 
minosity and a, though good for the Lyc flux of 
the regions on the envelope, should not hold for 
Ha since Ha does not measure the total emit- 
ted flux in this case. A first-order estimate of 
the appropriate index relating both the total Ha 
emitted flux and the Lyc flux can be obtained 
from the scaling model in Section 3.1. If the Ha 
luminosity of an H 11 region above the Stromgren 
transition is Ln^, the Lyc continuum luminosity 
(in units of the transition luminosity values) is: 

Lhyc = -^Ha"^ (22) 



and incorporating this relation into the virial pre- 
diction, eq. (9), we have 

iHa=pa4/^"=pa", (23) 

where p is a constant, and n = 4/ea. We 
can see that the condition for density bounding, 
e a > 1 implies that n < 4, and that the observed 
values for MlOO and MlOl: n ~ 2.5 are in agree- 
ment with this condition, yielding a value of a e 
close to 1.6. Equation (23) gives, in principle, a 
more reliable method for measuring a e than the 
LF slope ratio in eq. (21), but it requires a much 
higher investment in observation and reduction. 
If we can combine the two tests for a number of 
galaxies, this will give a useful quantitative check 
on the density-bounding hypothesis. 

3.3. The stability of the transition lumi- 
nosity 

In Table 1 we list the luminosities Lstr of the 
peaks in the observed LF at the transition. The 
mean value of Lstr is 4.05 x 10^^ erg s~i, which 
is 1.4 X 10^^^ photon s^^.This is the equivalent 
of ten main sequence stars of type 07, or 2 of 
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type 03 (Vacca, Garmany, & Shull 1996), which 
should be raised by a factor of 2.23 to aUow for 
the Lyc/Ha ratio. The rms scatter is surpris- 
ingly low: 6.2 % of Lsti, i-e., 0.07 stellar mag- 
nitudes. It is important to note here that if the 
peaks were induced by some effect of limited an- 
gular resolution, the rms scatter would be 0.4 
magnitudes, due to the range of distances of the 
galaxies observed, a powerful reason to reject this 
possibility, recently proposed by Plcuss & Heller 
(2000). It is interesting to compare this with 
the dispersion of the integrated luminosities of 
the same galaxies, taken from the Tully catalog 
(1987), which show a scatter of 0.15 mag. Al- 
though the sample is too small to draw a pow- 
erful conclusion, the promise of the glitch as a 
distance indicator is supported by the fact that 
its scatter is significantly lower than that of the 
absolute luminosities of the host galaxies in the 
sample under observation. 

Our sample is still too small to determine 

whether the transition can be used as a stan- 
dard candle on extragalactic scales, but it ap- 
pears a priori to satisfy the four criteria pro- 
posed by Aaronson Sz Mould (1986) for a cos- 
mological standard candle. It should: (a) ex- 
hibit small, quantifiable dispersion; (b) be mea- 
surable in enough galaxies to be calibrated lo- 
cally; (c) have a well-defined physical basis, and 
(d) be bright enough to be used well into the re- 
gion where the Hubble flow predominates. The 
observed peak in the Ha LF for H ii regions in 
disk galaxies promises to satisfy all four crite- 
ria. We have given here the empirical case for 
(a), (b), and (d), and very reasonable backing 
for (c). A criterion not brought out by Aaron- 
son and Mould is that the observable should 
be permanently available to observe. This is of 
course not satisfied by supernovae, and although 
permanency is not an absolute requirement, the 
Stromgren transition appears to be a high lumi- 
nosity and also a permanent feature of disk (and 
probably irregular) galaxies. The initial work 
presented here clearly requires support from a 



larger sample of galaxies. 

Although a purely empirically determined value 
of Lstr could, given its low scatter, allow it to be 
used clS 8j standard candle, it is clearly desirable 
to establish the underlying physics more firmly 
than we have been able to do here. Nevertheless, 
if we can assume the transition hypothesis main- 
tained in the present paper it does give a clear 
physical explanation for the low scatter in Lstr- 
Starting from eq. (10), and assuming a constant 
cloud density (which, as explained above in Sect. 
2.2 appears valid for clouds up to and includ- 
ing those with Lho = -^Str) we can rewrite the 
equation in terms of the mass, Mj, of the ionized 
cloud: 



Mi = k'M^, 



(24) 



where k' is a constant. The stellar mass, M*, 
and the mass, Mci, of the whole placental cloud 
are related by equation (11). At the transition, 
the mass of the ionized cloud just equals the 
whole cloud mass, and calling this mass Mstr we 
can find an expression for Mstr by combining (24) 
and (11): 



(25) 



where j is the constant of equation (11). 

If we take a = 1.5 = a e, the fractional varia- 
tion in the mass of a cloud for a given luminosity 
will be: 



(26) 



However, in any given galaxy, the transition 
peak in the luminosity function is formed by a 
significant number of regions, which we will spec- 
ify as N, so that the variance in the peak lumi- 
nosity will be given by: 



dMstr 
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Since in the galaxies observed N is of order 
50, N^^^"^ is of order 0.15, and it is this collec- 
tive property of Lgtr which makes it such a stable 
index as opposed to say the luminosity of the ra*^ 
brightest individual H ll region. We would ex- 
pect the main variable influencing both k' and 
j to be the metallicity. The parameter j de- 
pends on the IMF, and recent results show, to 
the surprise of some, that the metallicity appar- 
ently hardly plays a role in determining the clus- 
ter IMF slope. The IMFs of stars in the Galaxy 
(Massey, Johnson, &; de Gioia-Eastwood 1995) 
and in the LMC (Massey et al. 1995), with 
masses greater than 10 Mq are found to have 
slopes close to the Salpeter (1955) value of -2.35, 
in spite of their large metallicity difference. In 
a recent review article, Elmegreen (1997) argues 
that this is because of the essentially fractal na- 
ture of the fractionation process leading to star 
formation. Whatever the physical cause, this in- 
variant property works in favor of an invariant 
luminosity for the Stromgrcn transition. 

The parameter k' relating the number of Lyc 
photons to stellar mass is of course a function of 
metallicity, but over a range between 0.5 solar 
and 2 solar, calculations by Garcia Vargas et al. 
(1995) show that the flux from an O star will not 
vary by more than a few per cent, and the quan- 
titative effect of different metallicities on the net 
absorption within the H ii regions themselves are 
also small, if we are looking at hydrogen ioniza- 
tion (the situation would be more metal depen- 
dent for helium). 

In practice, to establish Lstr as a standard 
candle requires further steps: calibration locally 
via Cepheid distances using the HST key project 
galaxies and extension of the LF data base in 
nearby galaxies to obtain a reliable variance for 
Lstr- The potential range in z would be limited 
by the angular sizes of the H II regions, rather 
than by their luminosities. A practical ground- 
based limit is given by the distance of an H II 
region with Lhq = ^str which subtends 0.5 arc- 
sec, i.e., 2x10"^ km s~^. The use of HST would 



push this limit back considerably. Since the vari- 
ance appears to be low, this extra secondary dis- 
tance standard may be of use, e.g., in probing the 
depth-structure of neighboring galaxy clusters. 

4. Escaping Lyc photons, and the diffuse 
Ha emission 

The origin of the photons which ionize the 
diffuse interstellar medium relatively far from 
H II regions, which does not contain observable 
sources of ionization has not been definitively as- 
signed. The OB stars in H II regions are in some 
sense a natural source of this ionizing radiation, 
but doubts have arisen for two reasons. Firstly 
can enough radiation escape from the immedi- 
ate environments of the OB stars, i.e., from the 
H II regions themselves, and secondly can the 
mean free path of these photons be sufficiently 
long for them to ionize the diffuse medium over 
the full extent of a galactic disk? These doubts 
have led to at least one alternative model, no- 
tably than by Sciama (1990) who proposed the 
decay of massive neutrinos as a means of produc- 
ing in situ Lyc photons in the diffuse ISM. The 
hypothesis of density bounding for the most lu- 
minous regions can resolve the first of the two 
doubts (the second would have to await detailed 
models of the inhomogeneous ISM for its resolu- 
tion) . We can estimate the flux leaking out of the 
luminous H II region population by extrapolating 
the Ha LF of the H II regions in a galaxy in the 
range Lhq < -^^str to the range Lhq? > ^str, and 
subtracting off the measured LF in this upper 
range. The difference should approximate the 
ionizing flux equivalent which escapes from the 
density-bounded regions. If this is greater than 
the measured diffuse flux, we have a prima facie 
case in favor of the model as an explanation for 
the latter. 

We have used NGC 157 as an example to show 
whether this hypothesis is worth pursuing, and 
Figs. 7 and 8 illustrate the procedure for this 
galaxy. To take into account H II regions below 
our detectability limit, we extrapolated our LF 
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below this limit, using the LF measured for M31 
by Walterbos &: Braun (1992), which goes down 
to log Lhq = 35 (erg s~^) thanks to the prox- 
imity of M31. We scaled the measured LF for 
NGC 157 at log Lho = 38 where our sample is 
statistically complete, and used the form of the 
LF from Walterbos & Braun to extrapolate to 
lower values. The comparison of the observed 
truly diffuse Ha flux with that available due to 
Lyc photons escaping from the density-bounded 
H II regions is: 

• Observed upper bound to the diffuse flux: 
2.3 X 10^^ erg s~^ (implied ionizing photon 
flux 1.7 X 10^^ Lyc photon s~^). 

• Computed escaping flux: 3.1 x 10^^ ergs~^ 
(2.3 X 10^3 photon s"^). 

This example shows that, on the density- 
bounding hypothesis, the flux escaping from the 
luminous H II regions is easily sufficient to yield 
the diffuse Ha observed in NGC 157, and in a re- 
lated study (Zurita, Rozas, & Beckman 2000) we 
are amplifying the sample treated in this way to 
see whether the model is as good at accounting 
for the diffuse Ha in a wider sample of galax- 
ies. Ferguson et al. (1996) showed that there is 
a good geometrical correlation between the po- 
sitions of H II regions in general and the diffuse 
flux, in the two galaxies they studied in detail. 
It will be interesting to see if this correlation is 
improved when the most luminous H II regions 
are used, as should be the case if our model is 
valid. At this stage, given the dumpiness of all 
H II regions, we would not postulate that Lyc 
photons cannot escape at all from regions in the 
range Lhq < -t'Str, only that the regions above 
Lstr should be contributing proportionally in ex- 
cess of their already high measured Ha output. 

It is interesting that the computed escaping 
flux in NGC 157 is in excess of the observed dif- 
fuse flux. This difference appears to be well out- 
side our estimates of random and systematic er- 
rors in all the quantities estimated, and suggests 



that the down-conversion of Lyc photons in the 
diffuse medium is by no means sufficient to soak 
up all the ionizing flux leaking out of the density- 
bounded regions. A good fraction of all the ioniz- 
ing photons emitted by the OB stars must escape 
completely from the plane of the galaxy. If this 
effect is true of the majority of disk galaxies, the 
implication for the ionization of the intracluster 
medium in clusters of galaxies is of considerable 
potential interest. 

5. Discussion 

We have used three types of observations to 
argue for a change in physical regime between 
ionization bounding and density bounding in H II 
regions at a characteristic Ha luminosity close to 
log = 38.6 erg s^^ which we have termed 
the Stromgren luminosity, Lstr. They are the 
change of slope, accompanied by a "glitch" in 
the Ha LFs of the complete H II region popula- 
tions in a set of observed disk galaxies, the quan- 
titative relation between internal velocity disper- 
sion, a, and Ha luminosity, Lhq, of the H ii re- 
gions on the virial envelope of the log LfjQ-log 
a distribution for the complete population in a 
galaxy, and the change in behavior of the Ha 
surface brightness of regions in the range ~ 
Lgtr. Of these, the first two are presented as ev- 
idence of the kind of transition proposed, while 
the third, though of interest, is not as precisely 
defined by the observations, and its cause is at 
this stage fully open to other explanations. In- 
deed in order to sustain our basic hypothesis we 
will clearly require spectroscopic observations, as 
specified below. However, the narrow range in 
the luminosity of the peaks of the LFs presented 
here means that this feature is of considerable 
interest on purely empirical grounds, clS Si poten- 
tial secondary standard candle. For this reason 
alone it is worth making a firm effort to under- 
stand the underlying physical behavior causing 
the peak and change of slope. 

In Sect. 2.1 we showed that two phenomena 
which could occur, and which at first thought 
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might be able to account for the LF observations 
(the effect of the overlapping of regions on scales 
commensurate with the luminosities around Lhq 
= i^str, and the effect of dust extinction within 
and enveloping the regions) will not give rise to 
the observed LF change, and are not in fact vi- 
able as explanations for the observations. We 
ought also to consider an effect brought out by 
MW, who examine the relationship between the 
properties of H ll regions within the Galaxy, and 
their OB stellar components. They present mod- 
els to explain the change in slope of the Ha LF 
at moderately high luminosities, which they at- 
tribute to an essentially statistical effect due to 
the discrete number of high-luminosity stars, and 
its relation with the total ionizing luminosity of 
an H II region. In their paper, MW discuss the ef- 
fects of envelopes of diffuse gas around giant H II 
regions, which are observed for the majority of re- 
gions, and which correspond to the "core-halo" 
structure noted in Kennicutt, Edgar, &; Hodge 
(1989). In the present study we have taken this 
structure as correct, and have not discussed it in 
further detail. The regions defined by our empir- 
ical limiting isophotal method contain the whole 
of the core and the halo, and are clearly bigger 
than regions defined only via their cores, or via 
their emission radii at radio wavelengths. It is 
important to note here that the change in the 
LF slope predicted by the MW models occurs at 
an Ha luminosity an order of magnitude less than 
that of the glitch we observe. Our observations 
here effectively disguise the MW break, because 
it occurs where our sample is incomplete, leading 
to the broad peak and decline to low luminosities 
which are statistical artifacts. Neither our mod- 
els nor our observations contradict the results of 
MW, but their work does not examine the LF 
gradient change at Lh^ = -^^Str- 

The changes which occur in the surface bright- 
ness and surface-brightness gradients of the re- 
gions in the luminosity range L ~ Lstr are clearly 
of interest physically, but it is not so easy for us 
to show that they are attributable to a transition 



to density bounding, although we do have signif- 
icant technical advantages in resolution and sig- 
nal to noise over the observers in the 1980s (for 
example the number of regions listed in Kenni- 
cutt, Edgar, k Hodge 1989 for NGC 7479 was 
67, while the in data yielding the internal bright- 
ness gradients used in the present paper we cat- 
alogued over 1000 regions in this galaxy (Rozas 
et al. 1999). Nor can the changes in parameters 
presented in Figs. 5a and 5b be readily attributed 
merely to angular resolution limitations, as sug- 
gested recently by Pleuss & Heller (2000), based 
on their study of MlOl. Their argument is based 
on a scaling of the results from that galaxy to 
a distance of 20 Mpc, but the distances of the 
galaxies in the present sample range from 7 to 
38 Mpc {Ho = 65 km s"^ Mpc"^), a factor of 
over 5 in distance. The observed increase in sur- 
face brightness and surface-brightness gradients 
in the range Lua > -^str cannot occur for regions 
with constant density and filling factor, as well 
pointed out by McCall (1999); density-bounded 
regions would not show surface-brightness gra- 
dients different from ionization-bounded regions 
if their mean densities and filling factors were 
invariant. Thus although we have used the ap- 
proximation of constant density and filling fac- 
tor in modeling the changes in properties across 
the transition, it cannot hold for regions at much 
higher luminosities than Lstr- The product of 
density and filling factor must rise in this range. 
We can examine a simple case in which two 
H II regions have the same Ha luminosity but 
one is ionization-bounded and the other density 
bounded; the factor relating the total ionizing 
luminosity of the density-bounded region to the 
fraction which is down-converted to Ha we term 
g. In a self-consistent scenario, g can represent 
the factor by which the product of the density 
and filling factor of the density-bounded region 
exceeds that in the ionization-bounded region. It 
is easy to show that the central surface bright- 
ness of the (assumed spherical) density-bounded 
region will be the greater by a factor 5^/^, and 
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the radial-brightness gradient, parametrized as 
the central brightness divided by the radius, will 
be greater by a factor g. This exercise does not 
prove that density-bounded regions will necessar- 
ily have higher surface brightnesses, and bright- 
ness gradients, only that the assumption of den- 
sity bounding is consistent with this condition, 
and hence consistent with our observations. In 
order to fulfill this requirement, either the den- 
sity, or the filling factor, or both, must be rising 
parameters above = -^^Str- While it is per- 
fectly plausible that more intense ionizing sources 
produce higher filling factors, and that the most 
massive clouds have higher mean densities, it is 
not obvious that these changes will not begin to 
occur below Lhq = i^strj it is not obviously 
a physical necessity that the surface-brightness 
increase occur just at the transition. The fact 
that it appears to occur at, or near, what we 
believe to be the transition cannot at this stage 
be taken as strong support for the hypothesis of 
density bounding. Direct, spatially resolved mea- 
surements of in situ electron densities via line 
intensity ratios will be required to further our 
understanding of this point. We must also be 
prepared to examine our data in the light of non- 
isotropic models for luminous H II regions such 
as the "chimney" hypothesis (Norman & Ikeuchi 
1989; Heiles 1990) discussed in more physical de- 
tail by Dove & Shull (1994), in which at a critical 
luminosity a region can break out of the denser 
disk gas into the halo. These models appear, 
however, to predict reduced surface-brightness 
gradients for regions observed in face-on galax- 
ies, and so are unlikely to offer an explanation of 
our observations. 

One of the most attractive aspects of our hy- 
pothesis is that it offers a promising scenario to 
account for the diffuse Ha from gas dispersed 
over the disk of a galaxy, outside the H ii re- 
gions. Two quite recent studies have pointed to 
the emission of ionizing flux from OB stars as a 
satisfactory explanation of the Ha emission from 
the warm ionized medium in disk galaxies. Oey 



& Kennicutt (1997) made a direct comparison of 
the rate of emission of ionizing photons from the 
OB stars in the Large Magellanic Cloud, with 
the rate of Ha emission from H ll regions, and 
concluded that up to 50 % of the flux emitted by 
the stars is not down-converted within the H II 
regions, and is thus available in principle to ion- 
ized the warm diffuse medium. They compare 
this with an estimate of 35 % of the radiation 
emitted by the diffiise medium in Ha and con- 
clude that the OB stars do put out sufficient flux 
in this galaxy. A study by Ferguson et al. (1996) 
of NGC 247 and NGC 7793 concluded that the 
integrated flux required to ionize the diffuse gas 
is some 40 % of the total Ha flux emitted by a 
galaxy, including its H ii regions. However, this 
calculation assumed, as have any such previous 
calculations, that the total flux of ionizing pho- 
tons from an H li region can be directly measured 
via its Ha flux. If the hypothesis presented in the 
present paper is correct, a significant fraction of 
the total ionizing flux liberated within the H II 
region population of a galaxy is not converted to 
Balmer radiation within the H II regions. This 
is true of the most luminous regions, with lu- 
minosities ranging up to 10^^ erg s~^, and the 
escaping flux can range up to rather more than 
10^° erg s"-*^ for the brightest regions. The inte- 
grated escaping flux for a galaxy can attain a few 
times 10^^ erg s~^, easily enough to account for 
the measured diffuse Ha flux, according to the 
observations of Ferguson et al. (1996) and of our 
own (Zurita, Rozas, k. Beckman 2000). Ferguson 
et al. (1996) pointed out that the emitting dif- 
fuse gas tends to surround H ii regions, and we 
can go further on this point, noting the geomet- 
rical correlation with the most luminous regions, 
which we will go further to quantify in a new 
study (Zurita, Rozas, & Beckman 2000). Finally 
if our hypothesis is valid, a significant fraction of 
all the ionizing photons emitted by the OB stars 
in a galaxy may escape completely from the disk 
into the halo and finally into the intergalactic 
medium. Two corollaries are of interest here: in- 
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tergalactic clouds, even relatively far from a ma- 
jor galaxy, may have their surfaee layers ionized 
by this escaping flux, and global star formation 
rates in disk galaxies, estimated from integrated 
Ha fluxes, may in fact be significantly larger than 
these estimates. Both of these points are worth 
pursuing theoretically and observationally. 

We should treat the density bounding-hypothesis 
for the luminous H II regions as speculative. Fun- 
damentally missing pieces of the puzzle are line- 
ratio tests, of the sort performed by McCall, Ryb- 
ski, & Shields (1985) comparing [O iii] to [O 
ll] emission-line intensities (in a density-bounded 
region, the [O ll] Stromgren sphere should under 
many conditions be larger in radius than the H lI 
region, so the [O in]/[0 ii] ratio should be sig- 
nificantly enhanced over its value in ionization- 
bounded regions. McCall, Rybski, &; Shields 
(1985) claimed one such detection, region(-606- 
1708) in NGC 598, but no others, and this de- 
tection has since been called into question (Mc- 
Call 1999). Technical advances of the type rep- 
resented by the TTF (Taurus Tuneable Filter, 
Bland-Hawthorn & Jones 1998) which permits 
imaging of complete galaxies in single emission 
lines with full redshift flexibility, will allow us to 
apply line-ratio tests over complete populations 
of fully imaged H il regions in galaxies such as 
the ones measured only in Ha for the present pa- 
per. Another prospect, but only for the nearest 
galaxies due to limitations on angular resolution, 
is the use of two-dimensional fiber-fed spectro- 
graphs to sample the emission across the full face 
of an H II region. Here we would trade complete 
spectral coverage for limitation on the number 
of H II regions. In both cases we would hope to 
enhance understanding of the internal physics of 
the gas within the H II regions, measuring as a 
function of position on a region its temperature 
and electron density, and applying tests such as 
the [O lll]/[0 ll] emission ratio to the question 
of density bounding. 

The apparent invariance in the luminosity of 
the "glitch" measured in the LF does offer a pos- 



sible refined secondary standard candle of high 
luminosity and constant presence for use well 
into the Hubble flow, independently of whether 
the interpretation we have placed on it in the 
present article is fully, or even partially, valid. To 
demonstrate its scope will require amplification 
of the number of objects studied to determine 
reliably, using local galaxies, the scatter in the 
observed feature, including any possible depen- 
dence on galaxy luminosity and type, and to cal- 
ibrate it using galaxies whose Cepheid distances 
have been determined in the HST key project. 

The William Herschel Telescope is operated 
on the island of La Palma by the Royal Green- 
wich Observatory in the Spanish Observatorio 
del Roque de los Muchachos of the Instituto de 
Astrofisica de Canarias. This work was partially 
supported by the Spanish DGICYT (Direccion 
General de Investigacion Cientifica y Tecnica) 
via Grants PB91-0525 and PB94-1107. This re- 
search has made use of the NASA/IPAC Ex- 
tragalactic Database (NED) which is operated 
by the Jet Propulsion Laboratory, California In- 
stitute of Technology, under contract with the 
National Aeronautics and Space Administration. 
Thanks are due to Andrew Cardwell for a series 
of computational tests related to the H ii region 
luminosity function. We are grateful to Chris 
McKee, and to the referee, Marshall McCall, for 
comments which have helped to improve this pa- 
per. 
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Fig. 8. — Ha image of NGC 157 illustrating how 
the regions of highest luminosity also dominate in 
surface brightness, and how the diffuse Ha in the 
disk is correlated spatially with these. The ma- 
jor fraction of the ionization of the diffuse warm 
medium in disks should be due to these extremely 
luminous, and also density-bounded regions. The 
external isophote is at an emission measure value 
of 7 pc cm^^. 
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Fig. 5. — (a): Central surface brightness, Sc, in Ha, vs. the cube root of the Ha luminosity, L^a^^^, 
for H II regions in four galaxies, selected for minimum spatial blending with neighboring regions. Below 

is a good linear fit to the data, predicted for ionization-bounded regions with 
a constant product of density, electron density and filling factor. Above this value there is a clear change 
of regime; Sc rises much more sharply with L^^^; the scatter also widens considerably. This change must 
be caused by a change in the interactive behavior gas-stars, and cannot be due only to changes in stellar 
content or evolution within the regions concerned, (b): Radial gradient parameter; the ratio of the central 
surface brightness of a region to its estimated radius, in linear units. Near constant gradients are found 
for log IjHa < 38.6, and strongly increasing values with increasing scatter, for higher luminosities. The 
curves show the ratios of observed to intrinsic gradient for the two extremes of seeing-limited resolution 
(0.8" and 1.5") during the observation. A set of reigns with intrinsically constant gradients would show 
observed gradients tending to a constant value at high Lho, rather than the steeply rising gradients 
observed. 




Fig. 6. — (a): Log Lhq vs. the velocity half- width, a, for the most intense emission components from 
the H II regions in MlOO, observed using the TAURUS Fabry-Perot spectrometer. The linear fit is to 
the upper envelope of the plot, which should represent the locus of virialized regions. The slope is not 
4, as predicted for ionization-boundcd models, but 2.6, in agreement with density-bounded models for 
regions with log Lhq > 38.6. The same slope was measured by Arsenault, Roy, &; Boulesteix (1990) using 
highest surface brightness rather than lowest cr as a selection criterion, in full accord with the hypothesis 
of density bounding, (b): Plot of log Lhq vs. log a for the principal emission components in MlOl, from 
a TAURUS data cube without absolute flux calibration. A similar envelope at high Lhq and low a is 
found as for MlOO, with a slope calculated at 2.55, in accord with the hypothesis that the most luminous 
H II regions in MlOl are density bounded. 
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